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ABSTRACT 

There are several observational campaigns under way to detect kilometer size foregromid Kuiper Belt 
Objects (KBOs) through their occupation of background stars. The interpretation of an occultation 
light curve, unfortunately, is affected by a geometric degeneracy such that one is unable to determine 
the KBO size independent of its distance. This degeneracy can be broken through a precise measurement 
of the relative velocity obtained from simultaneous observations of individual events. While an array of 
telescopes spread over an area of few square kilometers can be employed, it is unlikely that the relative 
velocity can be measured to the required accuracy to help break this geometric degeneracy. The presence 
of diffraction fringes in KBO occultation light curves, when projected sizes of occulted stars are smaller 
than the Fresnel scale, improves determination of size and distance significantly. In this regard, we 
discuss the potential role of a dedicated satellite mission for KBO occultation observations. If dwarf 
stars at the V-band magnitudes of 14th and fainter can be monitored at time intervals of 0.1 seconds 
with normalized flux errors at the level of 1%, the occultation observations will allow individual KBO 
sizes and distances to be determined at the level of a few percent or better. 

Subject headings: occultations — Kuiper Belt objects — minor planets — solar system: formation 



1. INTRODUCTION 

The Kuiper Belt Object (KBO) population may provide 
useful clues to understand the formation and evolution of 
our Solar system and many others. The presence of KBOs 
beyond the orbit of Neptune is now well established; cur- 
rent estimates suggest that there is a total of order 10^ 
objects with sizes greater than 100 km (see Luu & Jewitt 
2002 for a recent review) . While sensitive photometric sur- 
veys may be capable of detecting KBOs with sizes greater 
than several tens of kilometers, the size distribution of 
KBOs with sizes of order a kilometer and below is use- 
ful for comparison with numerical models (Kenyon 2002). 
While imaging surveys are unlikely to detect kilometer- 
sized KBOs, as the expected R-band magnitude is of or- 
der ~ 30 or fainter, occultation surveys provide a much 
needed probe to extract population statistics at this low 
end of the size distribution. Unlike significant dependence 
on distance associated with reflected light, occultation ob- 
servations (e.g., Elliot & Olkin 1996) can also easily extend 
statistics of minor bodies beyond the current range of dis- 
tances encountered. 

In previous discussions, the use of occultation observa- 
tions to understand the KBO population was considered 
(Bailey 1976; Brown & Webster 1997; Roques & Moncu- 
quet 2000; Cooray & Farmer 2003). Due to the potentially 
large number of KBOs at the low end of the size distribu- 
tion, the probability for occultation is at the level of ~ 0.1 
per star per year. This event rate can be compared to that 
of galactic microlensing surveys where the optical depth 
is at the level of 10~^ requiring monitoring observations 
of order million stars or more. The KBO occultations, 
however, have a significant disadvantage in that the time 
duration of events is typically at the level of a few tenths 
of a second or below. This puts severe constraints on the 
instrumentation side for reliable detections. 

In Cooray & Farmer (2003), we explored the use of 
small telescopes arrays, such as the one planned by the 



Taiwanese- American Occultation Survey (TAOS; Liang et 
al. 2002), and showed that these surveys can be used 
to constrain the size distribution of KBOs between few 
100 meters and 10 km. Such statistical constraints as- 
sume prior knowledge on the KBO distance and veloc- 
ity distributions. While this approach may be useful for 
first generation surveys, for precise knowledge on the KBO 
population, it is useful to constrain the KBO size distri- 
bution independently. A step towards this direction is to 
extract parameters, such as the distance and size, from 
individual light curves. Unfortunately, such an extraction 
is affected by a geometric degeneracy between the KBO 
size and its distance such that one cannot determine one 
of these parameters independent of the other. This degen- 
eracy, however, can be partly broken with a measurement 
of the relative velocity associated with the event and, in 
general, require observations with multiple telescopes for 
this purpose. 

The occultation light curves are expected to show effects 
related to diffraction when events involve background stars 
whose projected sizes at KBO distances are below the Fres- 
nel scale. In addition to enhancing the event rate (Roques 
& Moncuquet 2000; Cooray & Farmer 2003), through the 
presence of fringes, the diffraction effect has another ad- 
vantage in that one can break the degeneracy associated 
with KBO size and distance with no independent measure- 
ment of the relative velocity. In practice, such a parameter 
estimation require flux measurements of background stars 
at the level of a few percent or better. An interesting 
possibility for this purpose involves a dedicated satellite 
mission for KBO occultations. 

The discussion is organized as follows. In the next Sec- 
tion, we will summarize important aspects related to KBO 
occultations, such as the normalized light curve. In Sec- 
tion 3, we discuss practical considerations related to KBO 
occultations and consider requirements for an array of sub- 
meter class telescopes based on the importance of a rela- 
tive velocity measurement. In the same Section, we also 



2 



1.5 



1.0 



u- 0.5 I- 

■D 
<D 
N 

re 
E 



0.0 



-0.5 



CO 

E 



1.0 



0.5 



0.0 



-0.5 



-1.0 



^ditf 

dA,„/dR. 

dA,.,/dR 

dA„/dD 

dA,.ydv 




(a) 






V 












R.=0.25km 


' ' ' 

\/ 






1 1 1 ■ 1 ■ ■ ■ 1 


(b) 


- 




- 








R.=6.0km 


1 ... 1 ... 1 





0.0 



0.2 



0.4 0.6 
t (sec) 



0.8 



1.0 



Fig. 1. — The fractional dimming of background stellar inten- 
sity during an occultation event by a foreground KBO of 10 km 
in radius with a projected stellar radius of 0.25 km (a) and 6 km 
(b). Note that the full light curve is symmetric about the central 
position (taken to be at t=0 seconds) and we only show one-half 
here. The Fresnel scale is taken to be 0.7 km, consistent with the 
value at a distance of 40 AU and an observing wavelength of 500 
nm. In addition to the light curve (solid line), we also show the 
derivatives of the light curve in terms of the physical parameters re- 
lated to the occultation. As the star size is increased, the diffraction 
fringes disappear and the information contained within the deriva- 
tives decrease. This results in a substantial lowering of the ability to 
extract parameters from occultations that involve larger projected 
stars, which also happen to be the brighter ones, at KBO distances. 



discuss the ability to extract information from individual 
light curves obtained at high precision with a space-based 
monitoring mission for KBO occultations. 

2. KBO OCCULTATIONS 

In general, one can characterize an occultation of a back- 
ground star by a foreground KBO with two parameters 
involving the duration, AT, of the occultation and the 
fractional flux decrease, A, of the background star dur- 
ing the occultation. In the geometric description of an 
occultation, these two observables are related via AT = 
(H-VA) /w, where D is the distance to the foreground 
KBO, V is the relative velocity, and is the angular size of 
the background star. The latter is what one obtains from 
stellar physics based arguments, such as through colors 
or a spectrum, or from interferometric observations (van 
Belle 1999). Since A = R^/Rl, with KBO size, i?, and pro- 



jected star size, D4>^), with occultation light curves, 

the combination of R/v and D/v are well constrained when 
prior information related to the stellar angular size is avail- 
able. A precise measurement of the relative velocity during 
the occultation event is needed to estimate R and D sepa- 
rately, otherwise a small KBO near by or a large KBO far 
away can be equally used to interpret a given light curve. 
The relative velocity measurement requires observations 
of the same event at two separate locations that happen 
to lie within the same occultation path. Since occultation 
observations are mostly sensitive to kilometer size KBOs, 
this can be achieved with an array of telescopes where in- 
dividual elements are spread over an area of few square 
kilometers or more. For comparison, a similar degeneracy 
exists for microlensing light curves involving the lens mass 
and its distance and to help break this degeneracy, at least 
one additional observation is required from a satellite at 
Astronomical-Unit distances (e.g., Gould 1995). 

In the case of KBO occultations, another possibility 
to break the degeneracy between distance and size in- 
volves the presence of diffraction fringes in occultation 
light curves. At KBO distances, the diffraction effect be- 
comes important especially when the projected star size 
involved with an occu ltation i s smaller than the Fresnel 
scale given by i?F = At - 40 AU, i?F 0.7 

km when A = 500 nm and is typical of projected stellar 
sizes, especially for main sequence stars with magnitudes 
around 14. With diffraction, the normalized flux variation 
as a function of the separation between the foreground 
KBO and the background star, r, is (Roques et al. 1987): 

_ R, 

^diff(i?*,i?, |r|) 



/(|r-r'|,i?). 
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The functions Un{x,y) are Lommel functions: 
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Compared to the geometric description, the diffraction 
fringes separate R, D and v, and is no longer restricted 
to combinations of R/v and D/v. Additionally, i?* can 
also be estimated through associated smoothing though 
to improve the determination of other three parameters, 
prior knowledge on R^ is preferred. 

3. PRACTICAL ASPECTS 

3.1. Parameter Extraction 

When considering a potential search for KBOs based on 
occultation observations and to understand relative mer- 
its of a ground vs. a space based survey, we first consider 
the ability to extract physical properties of the KBO pop- 
ulation from individual occultation light curves. For this 
purpose, we make use of the Fisher information matrix: 

(4) 
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Fig. 2. — Constraints on the KBO radius and distance, assumed 
to be at 40 AU with a radius of 5 km, based on the occultation 
light curve. In (a), we consider a ground based survey with 0.1 
sec timing capabiUty and normaUzed flux errors of 10% at the 12th 
magnitude corresponding to a stellar projected size of 3 km. With no 
prior information on the stellar angular size, one cannot constraint 
any of the KBO parameters. As shown with parallel lines, with 
prior knowledge on the projected star size at an accuracy of 10%, 
one constraints a combination of distance and radius while the two 
are not independently determined. To break this degeneracy, prior 
information on the relative velocity, say from a telescope array, is 
needed and with velocity known to an accuracy of 5% (20 it 1 km 
sec~^) this leads to an useful determination of the radius but not 
the distance (middle ellipse). To further improve, we consider an 
optimistic case involving a 1% known value on the projected star 
size and 5% normalized flux errors (inner ellipse). In this case, the 
radius can be determined to a level of 10% while the distance still 
remains uncertain at a level of 25%. In (b), we consider a potential 
space based scenario with timing at 0.1 seconds and a precisely 
known wavelength for observations. Note the change in y-axis range 
between (a) and (b). With a star of projected size 0.1 km, known 
to 5% accuracy apriori, and with normalized flux errors at 5%, one 
breaks the degeneracy in KBO distance and size such that the radius 
is extracted at a level better than 2% and the distance is determined 
at a level of 5% (outer ellipse). Further improvements are obtained 
when the star size is known more precisely: middle ellipse shows the 
case with 1% while the inner ellipse shows an optimistic case with a 
precisely known value for the projected size of the background star 
and 1% error on normalized flux. 



whose inverse provides the optimistic covariance matrix for 
errors on the associated parameters (e.g., Tegmark et aL 
1997). In Eq. 4, L is the hkehhood of observing data set x, 
in our case the hght curve, given the parameters pi . . .pn 
that describe these data. Following the Cramer-Rao in- 
equality (Kendall & Stuart 1969), no unbiased method can 
measure the ith parameter with standard deviation less 
than (Fii)~^/^ if all other parameters are exactly known, 
and less than [(F^-'-)^^]^/^, from the inverse of the Fisher 
matrix, if other parameters are to be estimated from the 
data as well. 

For the present case involving KBO occultation obser- 
vations, we assume light curves will be sampled with time 
intervals of 5T (< AT). To describe a real observational 
situation, we also assume noise in the light curve with a 
variance, a^{Adis), consistent with the instrumental de- 
scription. We can now write the Fisher matrix associated 
with KBO occultation light curves as 

- a^A^ift) dp, dp, • 

Here, dA^ig/dpi denotes derivatives of the light curve, 
^diff (Eq. 1), with respect to the parameter pi. 

In Fig. 1, we show light curves and their derivatives 
with respect to the four parameters. In order to show 
the importance of diffraction in breaking degeneracies, we 
show these derivatives as a function of the projected stellar 
size. As this size is increased, diffractions fringes begin 
to disappear and the derivatives with respect to physical 
parameters become smaller such that the ability to extract 
parameters from the light curve is reduced. Similarly, with 
increasing projected size for the background star, the light 
curves resemble that of a geometric occultation and the 
discussion related to the presence of degeneracies apply. 

To show this explicitly, we now consider parameter ex- 
traction and illustrate in terms of the KBO distance and 
radius, by marginalizing over other two parameters involv- 
ing the projected star size and the relative velocity. In the 
case of bright star occultations, as expected from a ground 
based survey, there is a significant degeneracy such that 
one cannot establish the size independent of its distance. 
We illustrate this in Fig. 2(a) with the two outer lines. 
To break this geometric degeneracy, one requires prior in- 
formation related to the relative velocity during the event. 
In the case of ground based surveys, this information must 
come from a second observation of the same event along 
the shadow path as it moves across the observer plane. 
Assuming few kilometer baselines, the best one can hope 
to achieve from ground, especially for slow occultations 
towards quadrature, is at the level of a 1 km sec~^. Us- 
ing this value as a prior, one can determine the radius to 
an accuracy of 10%, while the distance still remains un- 
certain at a level of 50%. This can be understood due 
to the fact that the uncertainty in the projected star size 
degrades the distance determination. To further improve, 
one should have precise knowledge on the projected sizes 
of background stars. The most significant improvements 
are associated with smaller, and thus fainter, stars such 
that diffraction effects become significant. 

In Fig. 2(b), we illustrate such a possibility using the 
same KBO parameters but now with a projected size for 
the background star of 0.1 km. This size roughly corre- 
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sponds to peak of the main sequence stars with V-band 
magnitudes between 14 and 16 (for example, Fig. 1 in 
Cooray & Farmer 2003). Now, one breaks degeneracies sig- 
nificantly such that the radius is known to better than 2% 
and distance is known to 5%. Note that this was achieved 
with no knowledge on the relative velocity. The diffrac- 
tion fringes allow the relative velocity to be determined 
together with KBO size and distance. We have assumed 
a single wavelength observation here. In a real situation, 
a finite width over wavelength is expected and this will 
lead to an additional smoothing of fringes. Assuming a 
25% width in wavelength of observations around 500 nm, 
we determined that errors on KBO size and distance are 
increased by a factor of 1.5 to 2 from the case with a single 
wavelength. This is under the assumption that the filter 
shape is known precisely. If the filter function is treated 
as an unknown, a few percent constraint on any of the 
physical parameters is not possible. 

3.2. Surveys 

In prior discussions (Roques & Moncuquet 2000; Cooray 
& Farmer 2003), the consideration was given to the ability 
of KBO occultation surveys to constrain the KBO size dis- 
tribution through statistics related to time duration and 
normalized flux change. Instead of this statistical ap- 
proach, we have now considered the possibility of using 
individual light curves to extract physical parameters of 
the KBO population. Unfortunately, as discussed, there 
are significant limitations on the ability of ground based 
observations to measure KBO size independent of its dis- 
tance. To break this geometric degeneracy, an array of 
telescopes, with capability to perform relative timing be- 
tween individual elements, can be used for simultaneous 
observations and to measure the relative velocity associ- 
ated with detected events. Since one expects KBOs with 
sizes at the level of few kilometers to be detected from 
ground based arrays, individual telescopes should be sep- 
arated at the same size. Though we have not considered in 
detail, there is an interesting optimization problem related 
to such an array: given the fixed number of telescopes, one 
can design the layout that will return the most informa- 
tion related to the KBO population. A simple estimate 
suggests circular arrays with separations both smaller and 
larger than the typical KBO size one expects to detect; the 
latter is to reduce systematics, as in the approach related 
to TAOS, while the former allows a relative velocity mea- 
surement. Due to small separations of order few hundreds 
of meters or below, it is unlikely that TAOS will be able to 
make reliable measurements of relative velocity associated 
with its occultation detections. For such surveys, infor- 
mation related to the KBO population can be extracted 
through statistics related to two observables involving the 
event duration and the normalized flux variation. 

Beyond TAOS, it will be useful to have an optimized ar- 
ray of small telescopes for occultation observations, though 
the information on the KBO population one extracts is 
still limited. We demonstrated this based on a calcula- 
tion related to the parameter extraction and showed that 
while KBO radii may be determined at the level of 10% or 
better, KBO distances will remain significantly uncertain. 
The significant improvements come with observations of 
occultations by KBOs of fainter stars, V-band magnitudes 
fainter than 14, such that diffraction effects are significant. 



The advantage here is that no separate relative velocity 
measurement is needed. For a successful ground based 
survey that attempts to make use of diffraction, observa- 
tions of fainter stars must be achieved with relative flux 
measurements down to few percent with timing capabili- 
ties around 0.1 seconds. 

To further improve, a dedicated space based satellite 
mission is helpful with the ability to monitor 14th mag- 
nitude or fainter stars with relative flux errors below 1% 
in integration times of order 0.1 seconds or better. The 
planned Kepler mission to detect transits of extra-solar 
planets across their central stars provide a basis to com- 
pare what will be required for KBO occultation observa- 
tions. The Kepler mission allows the detection of rela- 
tive fiux variations at the level of 10~^ on time scales of 
6 hours when 12th magnitude dwarf stars are involved. 
Scaling their expectations, we find that in an integration 
time of 0.1 seconds, at the 14th magnitude, the photon 
noise will lead to Poisson noise contributions to relative 
flux measurements of order 1%, which is the level required 
for a KBO survey. Since the probability for occultations 
is relatively high, the data rate can be controlled by a pri- 
ori selecting few hundred to thousand stars for monitoring 
purposes. 

To conclude, considering all possibilities, we suggest 
that a dedicated mission for KBO occultation observations 
may provide the best option to understand the KBO pop- 
ulation at kilometer sizes and below. Unlike observations 
of refiected light, occultations allow detections of KBOs, 
in principle, independent of its distance. The expected 
large number of KBOs at the low end of its size distri- 
bution leads to a significantly high occultation event rate 
and motivates ground-based attempts for KBO searches. 
While light curves are likely to be affected by a geometric 
degeneracy, diffraction provides a useful approach to break 
this degeneracy and improve KBO size and distance mea- 
surements. The planned TAOS survey provides an initial 
step in KBO searches while significant improvements can 
be expected with a well planned space-based mission. 
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